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. Abstract 
> ■ 

■ We analyze renormalization and the high temperature expansion of the 

■ one-loop effective action of the space-time noncommutative cp'^ theory by 
using the zeta function regularization in the imaginary time formalism (i.e., 
on S"^ X M^). Interestingly enough, there are no mixed (non-planar) contri- 
butions to the counterterms as well as to the power-law high temperature 

! asymptotics. We also study the Wick rotation and formulate assumptions 

' under which the real and imaginary time formalisms are equivalent. 



o 



K> ■ 

u 1 Introduction 



Since a number of review papers have been published recently (see [T]) it is 
not necessary to repeat here the motivations for studying noncommutative (NC) 
field theories. Most of the previous previous works, see e.g. |2l El IH E] and 
references therein, on finite temperature NC theories analyzed the case of space- 
space noncommutativity (with very few exceptions [6], [7]). Indeed, the case of 
space-time noncommutativity is most problematic because of the difficulties with 
unitarity and causality which were discovered some years ago [S], [HI HO] ■ These 
difficulties have not been completely resolved up to now. Space-time NC theories 
with compact dimensions exhibit an interesting phenomenon of discreteness of 
time [n]. 
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The main purpose of this paper is to develop certain aspects of the Euchdean 
space formahsm in space-time NC theories, including the renormalization, the 
transition from real to imaginary time, and the high temperature asymptotics. 

We shall start our work with analyzing the one-loop divergences in the Eu- 
clidean NC (j)^ on 5^ xM^ to make sure that the theory which will be discussed later 
does exist at least at the leading order of the loop expansion. We shall use the 
zeta-function regularization [121 [IS] and the heat kernel technique [T^ [T5l [T6] . In 
the context of an NC field theory the heat kernel expansion was first obtained for 
the operators which contain only left or only right Moyal multiplications p!71 ITS] . 
Such operators were, however, insufficient to deal with some physical applica- 
tions, like, for example, the 0*^ theory. The heat kernel expansion for generalized 
Laplacians containing both left and right Moyal multiplications was constructed 
in [19] on the Moyal plane and in [20] on the NC torus. Non-minimal operators 
were considered in [21j. We would also like to mention the calculations [22j of 
the heat kernel expansion in the NC (j)'^ model modified by an oscillator-type 
potential. 

To avoid unnecessary technical complications we shall study exclusively the 
case of pure space-time noncommutativity, i.e., we put to zero the NC param- 
eter with both indices in the spatial directions, ^■'^ = 0. We shall calculate 
the heat kernel coefficients with n < 4. It will appear that the coefficients 
02 and 04 look very similar to the commutative theory, but 03 is given by a 
complicated non-local expression. Fortunately, odd-numbered heat kernel coeffi- 
cients do not contribute to one-loop divergences at four dimensions in the zeta- 
function regularization. The model will turn out to be one-loop renormalizable 
with temperature-independent counterterms. 

Of course, we do not expect this model to be renormalizable at all loops. 
There are well-known problems related to the so-called UV/IR mixing [23] which 
should also be present in our case (though, maybe, in a relatively mild form since 
one of the NC directions is compact). To make the finite temperature NC 0^ 
renormalizable to all orders one should probably make it duality covariant [21] 
or use a bifermionic NC parameter [22] . 

An approach to finite temperature theories on static backgrounds based on 
the zeta-function regularization was developed long ago by Dowker and Kennedy 
|26j . In particular, they established relations between spectral functions of a 
3-dimensional operator which defines the spectrum of fiuctuations and the high 
temperature asymptotics of the free energy. In our case, due to the presence 
of the space-time noncommutativity, such 3-dimensional operator becomes fre- 
quency dependent even on static backgrounds. Therefore, eigenfrequencies of 
fiuctuations are defined by a sort of a non-linear spectral problem. Fortunately 
for us, a technique which allows to work with finite temperature characteristics of 
the theories leading to non-linear spectral problems has been developed relatively 
recently in the papers [27] • These papers were dealing with the thermodynam- 
ics of stationary but non-static space-times, but, after some modifications, the 
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approach of [27] can be made suitable for space-time noncommutative theories 
as well. By making use of these methods we shall construct the spectral density 
of states in the real-time formalism and express it through the heat kernel of a 
frequency dependent operator in three dimensions. Then, by using this spectral 
density, we shall demonstrate that the Wick rotation of the Euclidean free en- 
ergy gives the canonical free energy. To come to this conclusion we shall need 
two assumptions. First of all, we shall have to assume that the spectral density 
behaves "nicely" as a function of complex frequencies. Although this assumption 
is very hard to justify rigorously, we shall argue that the behavior of the spectral 
density must not be worse than in the commutative case, and we shall also sug- 
gest a consistency check based on the high temperature asymptotics. There is 
no canonical Hamiltonian in the space-time NC theories. Therefore, we have to 
assume that the eigenfrequncies of quantum fluctuations can replace one-particle 
energies in thermal distributions. This assumptions cannot be derived from the 
first principles of quantization basing on the present knowledge on the subject, 
but we can turn the problem around: the very fact that the Wick rotation of the 
Euclidean free energy leads to a thermal distribution over the eigenfrequences 
supports (a rather natural) guess that the eigenfrequences are the energies of 
one-particle excitations. Let us stress, that the calculations we shall perform in 
the Euchdean space do not depend on the assumptions described above. 

We shall also use the heat kernel methods to calculate the high temperature 
asymptotics of the Euclidean effective action assuming that the background field 
is static. As in the case of the counterterms, there are non-planar contributions. 
The asymptotic expansion does not depend on the NC parameter (provided it is 
non-zero) and looks very similar to the commutative case. 

This paper is organized as follows. In the next section we consider one-loop 
renormalization of NC 0^ on x M'^. Sec. |3] is devoted to the Wick rotation. 
High temperature asymptotics of the effective action are calculated in sec. ID 
Some concluding remarks are contained in sec. [51 

2 Noncommutative quantum field theory on x 
2.1 Basic definitions and notations 

Let us consider a scalar 0^ model on NC S"^ x M^. The scalar field is periodic with 
respect to the compact coordinate. We use the notations (x^^) = {x, x"^) = {x\ x^), 
where is a coordinate on S^, < t < f3. Similarly for the Fourier momenta 
we use k = (k, k^), k^ = I G Z. 
The action reads 

S = ll' dx' £^ d'x ((5,0)2 + m202 + ^ 0t) , (1) 
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where the (j)'l = (j)'k(f)'k(f)'k(j). Star denotes the Moyal product 

A * h{x) = exp [^e^'dpy) fi{x)f2{y)\y=.. (2) 

To simphfy the setup we assume that 9^^ = 0, but some of 6'^^ ^ 0, i.e. we have 
an Euchdean space-time noncommutativity. 

We wish to investigate quantum corrections to ([1]) by means of the background 
field method. To this end one has to spht the field into a classical background 
field if) and quantum fluctuations, = + b^f. The one-loop contribution to 
the effective action is defined by the part of ([T]) which is quadratic in quantum 
fluctuations: ^ 

S[ip,5ip] = - [ dx^ [ d^x 6<f{D + m^)6<f, (3) 
2 Jo Jr3 

where the operator D is of the form (cf. [281 [H]) • 

D = -{d,d^ + E), (4) 

with 

E = -f (L((^ ^^) + ^^) + L{ip)R{^)) . (5) 
o 

The one-loop effective action can be formally written as 

W = ^lndet{D + 771"^). (6) 

This equation still has to be regularized. To make use of the zeta-function regu- 
larization we have to define the heat kerne0 

= Tr(e-*^-e*^o) (7) 

and the zeta function 

C(s, D + m^) = Tr {{D + m^" - + m^)-') . (8) 

Here Tr is the L2 trace. In both cases we subtracted the parts corresponding to 
free fields with Dq = —d^j^d^ to avoid volume divergences. 
The regularized one-loop effective action is defined as 

Ws = --/i'^ -^e-'"^'K{t, D) = --/i^T(.)C(., D + m% (9) 

where s is a regularization parameter, is a constant of the dimension of mass 
introduced to keep proper dimension of the effective action. The regularization 

better name used in mathematics for this object is the heat trace, but here we use the 
terminology more common in physics. 
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is removed in the hmit s ^ 0. At s = the gamma-function has a pole, so that 
near s = 

Ws = ~(^--^E + In l/^ C(0, D + m^)~ ^('(0, D + m^), (10) 

where 7_b is the Euler constant. 

Let us assume that there is an asymptotic expansion of the heat kernel as 
t +0 

oo 

K{t,D) = Y,i^-~'^"^n{D). (11) 

n=l 

Such an expansion exists usually (but not always) in the commutative case. On 
NC S"^ X 'E? the existence of f llip will be demonstrated in sec. 12.21 For a Laplace 
type operator on a commutative manifold all odd-numbered heat kernel coeffi- 
cients a2fc-i vanish. (They are typical boundary effects). As we shall see below, 
on NC X the coefficient a-^ ^ 0. The coefficient oq disappears due to the 
subtraction of the free-space contribution in ([7]). 

The pole part of Wg can be now expressed through the heat kernel coefficients. 

C{<d,D + m^) = -m^a2{D) + ai{D). (12) 

Note, that odd-numbered heat kernel coefficients a2p-i{D) do note contribute to 
the divergences of Ws- 



2.2 Heat kernel expansion on x 

Let us consider the operator 

D = -{dl + E), E = L{h) + R{n) + L{k)R{r2) (13) 

on 5*^ X M^. This operator is slightly more general than the one in (jl]). The 
potential term is reproduced by the choice 

h = ri = -^(fi-kip , /2 = -r2 = wfv^. (14) 

D V 

We are interested in the asymptotics of the heat trace (171) as t ^ +0. To calculate 
the trace we, as usual, sandwich the operator between two normalized plane 
wave^ , and integrate over the momenta and over the manifold = S*^ x M^. 

K{t-D) = —^idk [ d4xe-*'="(e-*^-e*^")e*'^^ (15) 



^Although we are working with a real field, it is more convenient to use complex plane waves 
instead of real functions sin(A:a;) and cos(fcx). For a complex field we would have a coefficient 
of 1 instead of 1/2 on the right hand side of Since D with is real, this is the only 
difference. 
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where we introduced the notation 

^dk = J2 [ (16) 

with ^4 = 2im/ (3, n G Z. To evaluate the asymptotic expansion of f|T5l) at t — +0 
one has to extract the factor e~*^ . 

Kit^D) = j d'^fdke-^'^ 

x{e^^{t{{d-ikf + 2ik^{d^-ik^) + E))-l)u, (17) 

where we defined 

= e-**=^*Fe**=^ (18) 

for any operator F. Next one has to expand the exponential in (fT7|) in a power 
series in E and {d — ik). As we shall see below, only a finite number of terms 
in this expansion contribute to any finite order of t in the t +0 asymptotic 
expansion of the heat kernel. We push all [d — ik) to the right until they hit e*'^^ 
and disappear. 

K{t, D) = j d'xjdke-""' (tE + ^([9^, [d,, E]] + E^ + 2tk'^[d,, E]) 

4t^ \ 
-—k^k''[d„[d.,E]] + ...J^ . (19) 

We kept in this equation all the terms which may contribute to a„ with n < 4. In 
the commutative case all total derivatives as well as all term linear in k vanish. 
In the NC case this is less obvious because of the non-locality, so that we kept 
also such terms. The commutator of 9^ with ii^ is a multiplication operator, 
e.g., [d^,m] = L{d^l), [d^,L{l)R{r)] = L{d^l)R{r) + L{l)R{d^r). Therefore, 
eq. ( 1191) contains multiplication operators of two different sorts: the ones with 
only left or only right Moyal multiplications, and the ones containing products of 
left and rieht Moyal multiplications. The terms of different sorts will be treated 
different Ija. 

The terms with one type of the multiplications are easy. We shall call such 
terms planar borrowing the terminology from the approach based on Feynman 
diagrams. They can be evaluated in the same way as in ^TJ [TS] . Because of the 
identities 

' d^x{R{r))k= [ d^xr{x) , I d^x{L{l))u= I d^xl{x) (20) 



•^Formally R{r) ~ L(l)i?(r), but a constant function does not belong to C°°{S^ x M?) since 
it does not satisfy the fall-off conditions. Consequently, the two sorts of the term discussed 
above indeed lead to quite different asymptotics at t ^ +0. 
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only the E and E"^ terms in (fT9!) contribute. It remains then to evaluate the 
asymptotics of the integral 

dk e-^''' = (Ant)-^ + e.s.t., (21) 



/3(27r)3j 

where e.s.t. denotes exponentially small terms, to obtain 

^planarp^ = (^Tt)-' [ d^x{h + n), (22) 



Non-planar (mixed) contributions require considerably more work. The typi- 
cal term reads 

T(/, r) = I d'xjdke-'^' {mR{r))k (24) 

with some functions r(x) and l{x). For example, taking / = I2 and r = tr2 the 
expression (p^ reproduces the first term (tE) in (fT9|) . Let us expand r(x) and 
/(x) in the Fourier integrals 



^(^) = ^T7^(^f^^'%')e-'^ (25) 

Then 

(L(l)R(r))k = -J—.idq idq' f (g)/(g')e<'?+«')"e-t'=^(^-'?')-^('''-'=)^(^+'=) , (26) 
f3{2TTyj J 

where 

kAq = e^'^k^q, . (27) 
Next we substitute (J26i) in (JMl) and integrate over x and g' to obtain 

T{1, r) = .^l-^^jdkjdqe'''\-q) r{q)e-'''^^ . (28) 

In our case k A q = O^^k^qi — kiq^). 

Next we study the integral over k. The sum over k^ is treated with the help 
of the Poisson formula 

J]/(2™) = -^/ f{p)e-^-'dp. (29) 



We apply this formula to the sum 

J2 exp{-tkl - tO'^^hqj) , (30) 

fc4 

which corresponds to the choice 

/b)=-P^-;^--^J (31) 
in (l29l) . The sum fl30|) is transformed to (after changing the integration variable 



(5 . f°° 



The integral over kj is Gaussian and can be easily performed. We arrive at 
^ '""^ = (4^/^^''''^ V 4^ J 

where 

/i(g) = !{-q) f(g), = ^^^^^^J. (34) 

In eq. fl33l) one can still put |^^| = thus returning to the commutative case. 
The limit |^| — >^ does not commute however with taking the asymptotic t ^ 0. 
From now on we assume \6\ ^ 0. Obviously, all terms in the sum over are 
exponentially small as t — +0 except for q^ = 0. 

T{l,r) = ^/ d^q^^^W (^_(^!!^Z±M^ h{0,q)+e.s.t. (35) 
Let us define two projectors 

n;f = iri = 5'^ - (36) 

and split q into the parts which are parallel and perpendicular to O'^^: q\\ = 
Hyg, qs_ = U^q. Then d^q = dq\\£q^, and {O^^qj + pnf = (|^|g|| + pnf. The 
asymptotics of the integral over q\\ can be calculated by the saddle-point method. 
For each n there is one critical point of the integrand corresponding to q\\ = q^f* = 
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—(3n/\6\. We expand h{Q, gy, gj_) about these critical points and take the integral 
over q\\ to obtain 

T(/, r) = X: / (mo, ^ \ + ^/^"(o, g!" \ gx) + ...), 

(37) 

where prime denotes derivative with respect to q\\ . This completes the calculation 
of small t asymptotics for T{l,r). Since both /(x) and r{x) are supposed to be 
smooth, their Fourier components /(g) and f(g) fall off faster than any power 
at large momenta, and each term in the asymptotic expansion is given by a 
convergent sum and a convergent integral. 

The expression ( 1371) is already enough to calculate mixed (non-planar) contri- 
butions to the heat kernel expansions from the terms inside the brackets in f|T9l) 
which do not contain k. (We shall do this in a moment). Regarding the terms 
which do contain the momentum k, for our purposes it is enough to evaluate the 
power of t appearing in front of such terms. One can easily trace which modifi- 
cations appear in the calculations fl2^ - flTTI) due to the presence of a polynomial 
of kf^. The result is: (i) we still have an expansion in t^/^, (ii) the terms with k 
do not contribute to the heat kernel coefficients with n < 4. In other words, 
the only relevant mixed heat kernel coefficient is generated by the first term in 
the brackets in (HM . and it reads 



mixed 
"3 



(D) = E / ^'^^^~2(0, -gj"), -gx)r2(0, gj^^, g^), (38) 



where we substituted the fields appearing in E (see eq. ( fT3|) ). Note, that this 
expression is divergent in the commutative limit \9\ 0. The coefficient a™™'^ 
is highly non-local. The structure of (1381) . especially the sum over n, reminds 
us of the heat kernel coefficients on NC torus for a rational NC parameter |20j . 
In this latter case there is a simple geometric interpretation in terms of periodic 
projections [22]. No such interpretation is known for the present case of x M.^. 
However, some similarities can be found to the works [TT] discussing discretization 
of the coordinates which do not commute with a compact coordinate. 

2.3 Renormalization 

Let us return to our particular model ([T]). First we summarize the results of the 
previous subsection and re-express the heat kernel coefficients a„ = a^^^'^^'^ + a™™*^ 
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in terms of the background field (p by means of (fT^ : 



a,{D) = j dW, (39) 



,2 



aJD) = I d^xifit, (41) 

^ ^ 167r2 36 J ^ ' 

where tilde is used again to denote the Fourier components. 

Next we substitute (139!) - (14T!) in (ITOl) and (1T2|) to the pole part of the regu- 
larized effective action 

- / i'^ {^^,nV . . (42) 

This divergent part of the effective action can be cancelled by an infinite renor- 
malization of couplings in ([T]) 

Sm^ = - — - -, 5q = -. (43) 
There can be, of course, also some finite renormalization which we do not discuss 



here. Our main physical observation in this subsection the renormalization 
does not depend on the temperature 1//3. 

Here some more comments are in order. It is a very attractive feature of 
the zeta function regularization that the non-planar non-local coefficient 03 (D) 
does not affect the counterterms. This coefficient will, however, contribute at 
some other places, like the large mass expansion of the one-loop effective action 
(see, e.g., [IS])- Moreover, a^i^D) can lead to troubles in different regularization 
schemes. For example, if one uses the proper-time cut-off at some scale A defining 
the regularized effective action by 

1 df 

Wa = -- -e-'^'K{t,D), (44) 

the coefficient 03 generates a linear divergence oc Aas{D), which has no classical 
counterpart and cannot be renormalized away in the standard approach. There 
is a subtraction scheme (that was used in quantum field theory on curved back- 
ground [30] and in Casimir energy calculations [31]) which prescribes to subtract 
all contributions from several leading heat kernel coefficients, including as{D) in 
four dimensions. In the case of two-dimensional scalar theories this heat kernel 
subtraction scheme is equivalent to usual renormalization with the "no-tadpole" 
normalization condition In the present case the heat kernel subtraction is, 
obviously, not equivalent to the charge and mass renormalizations given by (H3|l . 
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We restricted ourselves to the case of pure space-time noncommutativity 9'^^ = 
0. However, one can try to make an educated guess on what happens for a 
generic non-degenerate 6"^*^. By comparing the heat kernel expansion obtained 
above with that on NC torus [20j and on NC plane with a non-degenerate 6'^'^ 
[T9] we can derive (rather unrigorously) the following rule: the presence of a 
non-compact NC dimension increases the number of the first non-trivial field- 
dependent non-planar (mixed) heat kernel coefficient by one as compared to the 
first non-trivial field-dependent coefficient in the commutative case. Indeed, in 
the commutative case the first such coefficient is 02- On the NC torus pO] (no 
non-compact dimensions) the first field-dependent mixed heat kernel coefficient is 
also 02- For the geometry studied in this paper (one non-compact NC dimension) 
this is 03. On an n-dimensional NC plane with a non-degenerate Of^" the first 
coefficient of interest is a„+2 [H]. We can expect therefore that the first mixed 
coefficient on x with a non-degenerate O^^" (three compact NC dimensions) 
will be 05. Such coefficient does not contribute to one loop divergence neither 
in the zeta-function regularization, nor in the proper-time cut-off scheme. Thus 
the situation in the generic case may be expected to better that in the case of 
a degenerate 6'^'^ discussed above. A similar conclusion has been made for the 
Moyal plane in p8] . 

As we have already mentioned above, the counterterms do not depend on 
the temperature. However, if one does the calculations directly on the zero- 
temperature manifold M^, there appears problems for a degenerate NC parameter 
|28] . Perhaps, compactification of one of the NC directions is a proper way to 
regularize these problems away. 

3 From imaginary to real time formalism 

The methods which allow to make correspondence between imaginary and real 
time formalisms in the case of frequency- dependent Hamiltonians were suggested 
in [27] and developed further in [53]. Here we briefiy outline these methods and 
discuss the peculiarities of their application to noncommutative theories. From 
now on we work with static background fields, doip = d^cp = 0. 

3.1 Spectral density in the real time formalism 

Let us consider a Minkowski space counterpart of the action ([1]). Our rules for 
the continuation between Euclidean and Minkowski signatures read ^4 — >■ tdo 
and 6^'^ —t6^^, where 6^^ is real, and 6^^do corresponds to O^^d^. We have, 
therefore, a real NC parameter in the Moyal product on both Euclidean and 
Minkowski spaces. These rules were applied, e.g., in [H], and they follow also 
from the requirement of refiection positivity [31]. As we shall see below, these 
rules also ensure consistency between the expressions for the free energy defined 
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in imaginary and real time formalisms. 

The wave equation for quantum fluctuations ipi^) over a static background 
reads 

(-^0^ + d]~m'- + Ri^') + Li^)Ri^))) ^(x) = 0. (45) 



The wave operator in fl45l) commutes with do. Consequently, one can look for the 
solutions ipi^ whose time dependence is described by ipcu{x) ~ e*'^^". They satisfy 
the equation 

iP{uj) + m^)i,^ = uj^^ , (46) 



where 
and 



P(cu) = + V{cu), Vico) = ^{^l + + ^+^„) (47) 



(4J 



Here we used the fact that left (right) Moyal multiplication of a function of 
by exp{iux^) is equivalent to a shift of the argument. 

From now on we consider the case of positive coupling g only. Then the 
potential V{u) is non-negative, V = {g/12){(p1_ + (p'^_ + (</?+ + (p-)"^) > 0. 

To define spectral density for the equation (H6|) we follow the works [27] . 
Consider an auxiliary eigenvalue problem, 

(P(A)+m2)^,,A = i^V.,A. (49) 



Obviously, the functions iIjuj,uj solve the equation (14611 . 

Our next step differs from that in [2?]. Let us restrict A to A < Aq for 
some Ao and put the system in a three dimensional box with periodic boundary 
conditions. Let us suppose that the size of the box is ^ 6'Ao, so that 99+ and 
if- are localized far away from the boundaries. In this case, the spectrum of the 
regularized problem can be considered as an approximation to the spectrum of 
the initial problem for the whole range of u. Later we shall remove the box, and 
the restriction A < Aq will become irrelevant. In the box, the spectrum of u in 
(H9|) becomes discrete, but, for a sufficiently large box, the spacing is small. The 
eigenvalues ^'^(A) depend smoothly on A not greater than Aq, and we can define 
the density of states as 

which can be used to calculate spectral functions of -P(A), where tilde reminds us 
that we are working with a finite-volume problem. For example, 

Tr:(e-*(P(A)W) ^ r dni^^^,^. ^^^^ 
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Here Tr denotes the L2 trace in the box. The potential V is non-negative. Con- 
sequently, there are no eigenvalues below m. 

The eigenvalues cu^ of the initial problem ( H6ll in this discretized setting ap- 
pear when the line v"^ = intersects ^'^(A). We can define the density of the 
eigenfrequencies uj^ by the formula 



dn{u!) 



(52) 



N 



Next, we would like to relate this density to fl50l) . This can be done by calculating 
derivative of the arguments of the delta function taken for u = X = u. We obtain 



dn{(jj) dn{uj,uj) 



where 



dn{i', X) 
rf(z/2) 



d(cj2 



rf(u;2 



N 



d{u 



N) 



d{X^ 



5{u' - u 



Nj 



(53) 



(54) 



This density admits an interpretation in terms of the heat kernel 

,-t(P(A)+m2) 



Tr 



\ 1 dP{X)\ 
2X dX J 



diy 



(55) 



Next we remove the box. Most of the quantities discussed above are divergent 
in the infinite volume limit. In order to remove these divergences we subtract the 
spectral densities corresponding the the free operator PQ + m"^ with Pq = —dj (not 
to be confused with P{0)). Then we perform the infinite volume limit. The limits 
of subtracted densities dn{u)/du, dn^u, X)/dX and dh^v, X)/dX will be denoted by 
p(a;), p(i/. A) and q{i', A), respectively. The following relation holds in this limit: 



Trs e 



-t{P(\)+n?) 



sub 



duj p{ijj; A) e 



where Trs is the L2 trace on M3 and 



Tr 



-t{P(\)+m^) 



3 e 



We also have the relation 



sub 



Trg ( e-*(^(^)+™') - e-*(-^.'+™') 



t(P(A)+m2) 



sub 



Q{v-X)e-"' dv 



(56) 



(57) 



(5^ 
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which, together with (!56|) . yields 



g{uj; \) = p{u; X) + - dxp{cr; \) da . (59) 

Jm 

To derive this formula one has to integrate by parts. Vanishing of the boundary 
terms is established by using the same arguments as in [27]. An infinite volume 
counterpart of fl5^ reads 

p{uj) = g{iu;iu). (60) 

An independent calculation of the spectral densities is a very hard problem. 
We shall view the equation fl56|) as a definition of the subtracted spectral density 
p(z/; A) through the heat kernel (an explicit formula involves the inverse Laplace 
transform). The other spectral densities £»(z/; A) and p(a;) are then defined through 
(EHD and (EDD. 

Relations similar to ( l56l) . ( l58l) . ( l59i) and ( l60i) were originally obtained in [27] 
for a different class of frequency dependent operators and by a somewhat different 
method. 



3.2 Wick rotation 

In this section we show that the Wick rotation of the free energy F defined 
through the Euclidean effective action coincides with the canonical free energy 
Fc. The methods we use are borrowed from [2^, but there are some subtle points 
related to specific features of NC theories. By definition, 

W{(3)=(3{F{(3) + S), (61) 

where S is the energy of vacuum fluctuations. 

Our renormalization prescription (1431) ia equivalent to the (minimal) subtrac- 
tion of the pole term (1421) in (llOp . Therefore, the renormalized one-loop effective 
action reads 

W = 'lj;\s=o{~^'X{s,D + m')), (62) 

where p^ := /x^e""^^. On a static background one can separate the frequency sum 
from the L2(M^) trace and rewrite the zeta function as 

C(s,D + m2) ^ J2TTs{{Luf + m^ + P{LUi))-'-iujf + m''-d])-') 
I 

/"OO 

= J2 dupEiiy;uJi){ujf + iyY'- (63) 

uji = 2711/(3. The spectral density p^; is defined for the Euclidean space NC 
parameter 9^'^. It is related to the real-time spectral density by the formula 

Pe{v; uj\e^^) = p{iy; tuj\ - iO^^) (64) 
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according to the rules which we have discussed at the beginning of sec. 13.11 We 
have aheady mentioned that the Wick rotation leaves the combination uj9 and 
the potential V invariant. Therefore, both densities coincide as functions of their 
arguments v and uj. However, we shall keep the subscript E to avoid confusion, 
but shall drop 9 from the notations for the sake of brevity. Next we use the 
formula 

^/M = £^£cot(^)/(^)d;^ (65) 

with the contour C consisting of two parts, C+ running from + oo to ?e — cxd 
and C- running from —le — cxd to —le + cxd, to rewrite the frequency sum as an 
integral. Then, by using the symmetry of the integrand with respect to reflections 
of z we replace the integral over C by twice the integral over C+ alone. Finally, 
we apply the identity 



to arrive at the result 

C(s, D + m')= (3Co{s, D + m^) + Ct(s, D + m^), (67) 

where 



CO oo 



Co(s, D + m^) = ^j du J pe{u; z){v' + z'Y'dz , (6 

m 

oo 

Ct(s, -D + vr?) = — du (b dz 



m Ca 



^ln(l-e*^") 
dz 



Pe{u;z){u' + z')-^ .(69) 



In commutative theories [27], the function (^y, which vanishes at zero temperature, 
represents the purely thermal part, while Co is responsible for the vacuum energy. 
In space-time NC theories there is no good definition of the canonical Hamiltonian 
and of the energy. Therefore, we have no other choice than to accept the same 
identities as in the commutative case, namely 

W) = -^^j-l=o~^'XT{s,D + m'), (70) 

£ = -\^^\s=ofi''Us,D + m''). (71) 

Actually, the definition of £^ is a rather natural one since it coincides with the 
renormalized Euclidean one-loop effective action on M''. However, as we have 
already mentioned in sec. l2.3l the renormalization in NC theories depends crucially 
on the number of compact dimension. Therefore, if one does the renormalization 
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directly in M^, one may need the counterterms which differ from (H3|) obtained 
on X R^. 

From now on we concentrate exclusively on Ft{I3) and Ct- We integrate by 
parts over z to obtain 



oo 

m C+ 



Ct{s) = / du d) dz ln(l — e 



dzpEji^] z) _ 2zspe{v] z) 



(72) 



To ensure the absence of the boundary terms we have to deform the contour C+ 
by moving its' ends up in the complex plane, so that e*'^^ provides the necessary 
damping of the integrand. We discuss the conditions on pE which make such 
deformations of the contour legitimate below. The integration by parts over u in 
the first term in the square brackets in (I72|) yields 



oo 

Ct(^) = ^ / / - ^'^') +^^2)s+i z) , (73) 



Til _ 

m C+ 



where 



Qe{v;z) = pE{y;z) -- I OzPe^cf] z)da . (74) 



The right hand side of (1731) is proportional to s. To estimate the derivative ds at 
s = in fITOl) one can put s = in the rest of the expression and use the Cauchy 
theorem after closing the contour in the upper part of the complex plane. The 
result is then given by the residue at z = iv. Next we make the Wick rotation 
of the NC parameter, so that pE^cr^ii') becomes p^a^u), and QEi^'^i^) becomes 
p(z/; z/) = p(z/) (cf. eqs. (15^ and (1^ ). Consequently, the Euclidean free energy 
is given by the equation 

oo 

F{f3) = ^ J dup{u)Hl-e-^n, (75) 

m 

which coincides with the canonical definition of the free energy Fc- 

The equality F = Fc is the main result of this section. To derive it we in- 
tegrated by parts and deformed the contour (7+. The integration by part over 
u is a safe operation, since for any fixed z the spectral density pe{i^, z) corre- 
sponds to the Laplace operator with a smooth potential. The absence of the 
boundary terms can be then demonstrated by standard arguments [27J based on 
the heat kernel expansion. The deformations of the contour are more tricky. To 
justify this procedure and application of the Cauchy theorem one has to assume 
that pEi^iz) can be analytically continued to the upper half-plane as an entire 
function of z. A rigorous proof of this assumption is hardly possible even in 
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more tractable cases of stationary commutative space-times [27]. We may argue, 
however, that this assumption is plausible. Consider pure imaginary values of 
z = IK. All deformations of the contour are done before the Wick rotation of the 
NC parameter 9. Therefore, ip± becomes complex, and = ^*_- The potential 
V{iK) remains real and positive. The background field (p is assumed to fall off 
faster than any power of the coordinates in real directions to ensure the existence 
of the heat kernel expansion. Such fields typically grow in imaginary directions 
(one can consider ip ~ e~'^^ as an example). Large positive potentials tend to 
diminish the spectral density thus preventing it from the blow-up behavior. It 
seems therefore, that the spectral density in our case should not behave worse 
than the spectral density in the commutative case. Another argument in favor 
of our assumption will be given at the end of the next section. 

The free energy (!75l) is expressed through a thermal distribution over the 
eigenfrequencies. In the absence of a well-defined Hamiltonian it is not guaranteed 
that this is the same as a thermal distribution of one-particle energies. This is 
a known problem of space-time NC theories which is beyond the scope of this 
paper. 

4 High temperature asymptotics 

As in the previous section we rewrite the regularized one-loop effective action ([9]) 
on a static background in the form 

Ws = -\^J^''ns) J2 Trs {uj' + Piuj) + m^)-;, (76) 

where the sum over the Matsubara frequencies is separated from the trace over 
the L2 functions on R^. As usual, we subtracted the free space contributions 
corresponding to = in P(a;) (which is indicated by the subscript "sub" in 
(176|) ). We remind that u = 27il/P, I E Z. To evaluate the high temperature 
(small P) asymptotics of Ws we split the sum in ( 17^ in two parts, 

Ws = W^=° + W^/°, (77) 

which will be treated separately. 
We start with which reads 

Wi=' = -lf^''T{s) Trs (P(0) + m^)-; = -^^^T(.) C(., P(0) + m'). (78) 

(The subtraction of free space contributions is included in our definition of the 
zeta function, cf. (IHl)). For each given uj the operator P{uj) is a three-dimensional 
Laplace operator with a scalar potential. All effects of the noncommutativity are 
encoded in the form of this potential. Therefore, as for all Laplace type operators 
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on R^, the zeta function in (!78l) vanishes at s = making ° finite. We can 
immediately take the hmit s — ^> to obtain the renormahzed expression 

W'='' = -^CiP{0) + m^). (79) 

In the rest of the frequency sum we first use an integral representation for the 
zeta function 



sub 



and then use a trick similar to the one employed in the previous section. Namely, 
we replace the operator in the exponential on the right hand side of flSUl) by 
a;^ + + P(A), expand each of the terms under the frequency sum in asymptotic 
series at c<j ^ oo keeping A fixed, and then put \ = uj. The result of this procedure 
reads 

a;^0 n=2 "'0 

= -^/i^^$^f;|c.|^— 2T(!^ + .)a„(P(a;) + m2).,b. (81) 

a;^0 n=2 ^ ^ 

Some comments are in order. Here we used again the fact that u;^ + m^ + P(A) for 
a fixed UJ is just a usual Laplace type operator in three dimensions. The large uj 
expansion of the heat trace in flSU]) is therefore standard and, as well as the usual 
large mass expansion is defined by the heat kernel coefficients (see, e.g., [IS]). On 
a manifold without boundary an asymptotic expansion (fTTj) with the replacement 
D ^ P + m? exists, and only even numbers n appear. The coefficient oq vanishes 
due to the subtraction, so that the sum in (ISTj) starts with n = 2. 

Now, we have to study the behavior of a„(P(cij) + m^)sub at large uj. These 
heat kernel coefficients are integrals over of polynomials constructed from the 
potential V{uj) and its derivatives. We can present them as 

an{P{uj) + m'U^ = an{P + m^)^^^ + a„(P(^) + m^)™-^ (82) 

where the first (planar) contribution contains all terms which are the products of 
either ip+ and its derivatives only, or of if - and its derivatives only (but not the 
products of (p+ and (p-). The rest is collected in the second (mixed) contribution. 
Obviously, no subtraction for the mixed heat kernel coefficient is needed. Due to 
the translation invariance of the integral over M^, the planar coefficient does not 
depend on uj. E.g., J d^xip"^ = J d^xip"^. Therefore, we drop uj from the notation. 

First, let us consider the mixed contributions to flSTl) . We assumed that 
the background field (f belongs to C^{S^ x M^). Therefore, it should vanish 



18 



exponentially fast at large distances. Since each term in a„(P(co') + TT^^^^ixed 
contains a product of at least one with at least one it should be of 
order 6*26"'"^''^^' for large u, where C2 and Ci are some constants. Ci is positive 
and characterizes the fall-off of if at large distances. C2 depends on n, on the 
amplitude of and on the functional form of a„. Up to an inessential overall 
constant the contribution of a mixed coefficient to flHT]) can be estimated as 

~ J2 k^-^e-^^l-^l ~ f^p-'-'l'-'^exp f_^^l£M\ (83) 

uij^O 1=1 VP/ 

(this sum is obviously convergent, so that one can remove the regularization 
parameter). If jS is small enough, namely (3 <^ Ci\6\, all terms in the sum fl551) 
are strongly suppressed, and the value of the sum can be well approximated by 
the first term 

^^«-3expf-^!^V (84) 



We conclude that the contributions of the mixed terms are exponentially small 
and can be neg lectecfl in the high temperature expansion of the effective action. 

Since the planar heat kernel coefficients do not depend on u we are ready to 
evaluate their contribution to (ISTi) by using precisely the same procedure as in 
Dowker and Kennedy [26] . 

oooo/„\ / o \ n-3+2s 



n=2 i=l V / V / 

= E r ( + ^ j C/^(2. + n - 3)a„(P + m^)^-^ (^|- j (85) 

with (ji being the Riemann zeta function. We remind that the index n in ( !85|) is 
even. The only divergence in ( !85l) is a pole in for n = 4. The corresponding 
term near s = reads 



/ 7") I 2 \ planar 

-a4(P + 



' /sub (4^^1/2 



1 



(86) 



On static backgrounds there is a useful formula which relates planar heat kernel 
coefficients of D and P, 

^niD + m'tr' = J^2^niP + m^)^— . (87) 

This formula follows from the analysis of the planar heat kernel coefficients pre- 
sented in sec. 12.21 and general formulae for the heat kernel expansion of Laplace 



*If one imposes a stronger restriction on the background requiring that (p is of compact 
support, then the mixed terms vanish identicaUy above certain temperature. 
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type operators [HI [151 IIHl • The coefficient P appears due to the integration of 
a constant function over the Euchdean time, and (47r)^/^ comes from different 
prefactors in the heat kernel coefficients in 3 and 4 dimensions. In particular, 
a,iD + m^)f2''' = -m'a2iD) + a,{D) = [3 / {AT^f'^a^{P + m')f2''' (let us re- 
mind that mixed 02 (-D) and ai{D) vanish). From (JTOl) and ( fT2i) we see that the 
divergence in the Euclidean effective is reproduced. This divergence is then re- 
moved by the renormalization of couplings fH3|) . After the renormalization, we 
collect all contributions to the effective action to obtain our final result for the 
high temperature expansion of the renormalized effective action 



W 



3/2 



3/3 



«2(P + m2)^— --C'(P(0)+m^) 



2 \ planar 



--a4(P + m;,^b 



7b + 2 In 



An 



(88) 



n=6 



n— 3 



It is instructive to compare the expansion ( !88ll to the one in the commutative 
case obtained by Dowker and Kennedy [26] (note, that the normalization of the 
heat kernel coefficients used in that paper differs from ours). We see that the 
C term is the same in both cases. The terms proportional to the heat kernel 
coefficients for the commutative case can be obtained from the expansion above 



by means of the replacement a„(P + m 



2 \ planar 
/sub 



,(P(0) + m2),,b. (In both 



cases subtraction of the free space contribution means simply deleting the highest 
power of m in standard analytical expressions jl6j). Let us write down explicit 
expressions for a couple of leading heat kernel coefficients. In the NC case we 
have 



2 ^ planar 
sub 



(47r)3/2 

1 



(47r)3/2 

The coefficients appearing in the commutative case are 



9^ 4 , 9 22 
—cp +-m if 

3o 3 



^9) 



a2(P(0) + m2 



a4(P(0) + 



1 



'sub 



'sub 



(47r)3/2 

1 

(47r)3/2 



j'i 9 2 

d x-ip , 



d^x 



9 A , 9 22 



(90) 



In both cases the corresponding heat kernel coefficients differ only by numerical 
prefactors in front of the same powers of if. 

The high temperature expansion does not depend on 6. In the limit 6^0 
(which is a trivial operation) one does not reproduce the corresponding expansion 
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in the commutative case. The hmits /3 — > and 6—^0 are not interchangeable 
because of the condition f3 <^ Cil^l which was imposed when studying the mixed 
contributions to the asymptotic expansion. 

In space-space NC theories a drastic reduction of the degrees of freedom in 
the non-planar sector above certain temperature was observed in [3]. This may 
be related in some way to the absence of non-planar contributions to the high 
temperature power law asymptotics in space-time NC theories found above. 

One can calculate the high-temperature asymptotics also in the real-time 
formalism. The key observation that the non-planar sector does not contribute 
remains valid also in this formalism. The heat kernel expansion in the planar 
sector has the standard form (though the values of the heat kernel coefficients 
differ from the commutative case). Therefore, one can repeat step by step the 
calculations of [27] and obtain an expansion for the free energy which is consistent 
with the expansion for the effective action derived above in the imaginary-time 
formalism. This is another argument in favor of the assumptions made in sec. 

m 

5 Conclusions 

In this paper we considered some basic features of the finite-temperature NC 
(j)^ theory in the imaginary-time formalism. We restricted ourselves to the case 
of pure space-time noncommutativity, 6^^ = 0. We used the zeta function reg- 
ularization and the heat kernel methods. Although we found highly non-local 
non-planar heat kernel coefficients, such coefficients do not contribute neither to 
the one-loop divergences, nor to the high temperature asymptotics. The theory 
can be renormalized at one loop by making charge and mass renormalizations, 
as usual. The counterterms do not depend on the temperature (as long as it is 
non-zero). We expect that the renormalization of this theory at zero temperature 
proceeds differently. The high temperature expansion of the one-loop effective 
action looks similar to the commutative case. The coefficients of this expansion 
do not depend on the NC parameter 6, but again, one has to assume that this 
parameter is non-zero. 

We have also studied relations between the imaginary and real time formu- 
lations. We found that the Wick rotation of the Euclidean free energy gives the 
canonical free energy modulo two assumptions. One assumptions about the be- 
havior of the spectral density on the complex plane is of technical nature. Another 
one is more fundamental, it concerns the interpretation of the eigenfrequnecies 
of perturbations as one-particle energies. 

An extension of our results to more general models containing gauge fields 
and spinors can be done rather straightforwardly. Gauge fields are particularly 
important to make connections to other approaches [HI Ej. Curved space-times 
will probably be difficult because of the problems with the heat kernel expansion. 
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Even in the case of a two-dimensional NC space with a non-trivial metric the 
heat kernel coefficients for a (rather simple) operator are known as power series 
in the conformal factor only [35]. 
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